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Driven vortices in confined geometry: the Corbino disk 
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The fabrication of artificial pinning structures allows a new generation of experiments which can probe the 
properties of vortex arrays by forcing them to flow in confined geometries. We discuss the theoretical analysis of 
such experiments in both flux liquids and flux solids, focusing on the Corbino disk geometry. In the liquid, these 
experiments can probe the critical behavior near a continuous liquid-glass transition. In the solid, they probe 
directly the onset of plasticity. 



1. INTRODUCTION 

In the mixed state of type-II superconductors 
the magnetic field is concentrated in an array 
of flexible flux bundles that, much like ordinary 
matter, can form crystalline, liquid and glassy 
phases. [^]j2j In clean systems the vortex solid 
melts into a flux liquid via a flrst order phase 
transition. jl| If the barriers to vortex line cross- 
ing are high, a rapidly cooled vortex liquid can 
bypass the crystal phase and get trapped in a 
metastable polymer-like glass phase. Q The di- 
versity of vortex structures is further increased 
by pinning from material disorder, which leads to 
a variety of novel glasses. Disorder-driven glass 
transitions are continuous, with diverging corre- 
lation lengths and universal critical behavior. |^,|| 

Of particular interest is the dynamics of the 
vortex array in the various phases and in the prox- 
imity of a phase transition. In the liquid phase 
the vortex array flows yielding a linear resistiv- 
ity. In the presence of large scale spatial inho- 
mogeneities, the liquid flow can be highly non- 
local due to interactions and entanglement. 
The correlation length controlling the nonlocal- 
ity of the flow grows with the liquid shear viscos- 
ity, which becomes large as the liquid freezes. At 
a continuous liquid-glass transition this correla- 
tion length diverges with a universal critical ex- 
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poncnt. In the solid phase the vortex array moves 
as a single elastic object under uniform drive, pro- 
vided the shear stresses are not too large. In the 
presence of strong spatial inhomogeneities, plas- 
tic flow occurs for large drives (or even for van- 
ishingly small drives in a glassy solid) and the 
response is always nonlinear. |^ The dynamical 
correlation length can be identifled with the sep- 
aration between free dislocations and diverges at 
a continuous melting transition. Probing spatial 
velocity correlations can therefore give informa- 
tion on vortex dynamics within a given phase, as 
well as on the nature of the phase transitions con- 
necting the various phases. 

As for ordinary matter, the shear rigidity of 
the vortex array can be probed by forcing the 
vortices to flow in conflned geometries. 
This type of experiments was pioneered by Kes 
and collaborators to study the shear rigidity of 
the two-dimensional vortex liquid near freezing 
in thin fllms. |llo|] More recently, patterned irra- 
diation of cuprate superconductors with heavy 
ions has made it possible to create samples with 
controlled distributions of damage tracks. [pl| We 
recently showed that an analysis of such exper- 
iments that combines an inhomogeneous scaling 
theory with the hydrodynamics of viscous flux liq- 
uids can be used to infer the critical behavior near 
a continuous glass transition, as well as to dis- 
tinguish between continuous transitions, such as 
that to a Bose glass, and nonequilibrium transi- 
tion to a polymer-like glass driven by interaction 
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Figure 1. The field profile, E{r){2'Kt/pfI), in the 
Hquid annulus of an irradiated Corbino disk. The 
inner and outer radii are i?i = 2iim and R2 = 
12/iTO, and ^ = 1/zm. The dashed hne is the 
1/r field profile in an uncorrelated liquid, with 
^ = 0. Inset: a sketch of the disk - the Bose glass 
contacts are not shown. 



and entanglement. ||] . Large scale spatial inhomo- 
geneities can also be introduced in the flow, even 
in the absence of pinning, by applying a driving 
force with controlled spatial gradients, as done 
recently by the Argonne group using the Corbino 
disk geometry. Q In this paper we illustrate the 
analysis of spatially inhomogeneous vortex mo- 
tion in both the liquid and the solid using the 
Corbino disk as a prototype of a novel class of 
experiments exploiting the effect of geometry to 
study the dynamics of vortex matter. 

2. LIQUID FLOW IN CHANNELS 

In the Corbino disk, with magnetic field along 
the disk axis (z direction), a uniform radial cur- 
rent density of magnitude J{r) = I / {2'Ktr) is in- 
troduced in the sample by injecting current at the 
center and removing it at the outer circumference 
of the disk (inset of Fig. 1). The current drives 
the vortices to move in circles about the axis. In 
the flux liquid, the dynamics on scales larger than 
the intervortex spacing, ao, is described by hy- 
drodynamic equations for the flow velocity v(r). 



which determines the local field from flux motion, 
E — nQ(j)QZ X v(r)/c, with no = l/ag. For simple 
geometries like the Corbino disk, where the cur- 
rent is spatially homogeneous in the z direction, 
hydrodynamics reduces to a single equation, j|J^ 

-7V-|-ryViv = ino0oz X J(r), (1) 
c 

where 7(T, H) is the friction, 77 (T, H) is the vis- 
cosity controlling the viscous drag from interac- 
tions and entanglement, and the term on the right 
hand side is the Lorentz force density driving flux 
motion. It is instructive to rewrite Eq. as an 
equation for the local field, 

-C'ViE + E = p;J, (2) 

with ^ = y^rf/j the viscous correlation length 
and pf = (no(/)o/c)^/7 the fiux flow resistivity. 
If the viscous force is negligible, Eq. (^) is sim- 
ply Ohm's law and the radial field is Eo{r) = 

(p///2^i)(l/0- 

To probe the viscous drag, it is necessary to 
force large scale spatial inhomogeneities in the 
flow. This may be achieved by suitable pinning 
boundaries. As an example, we imagine selec- 
tively irradiating a cylindrical central region and 
an outer annular region of the disk to obtain the 
structure sketched in the inset of Fig. 2. Here 
the vortices in the heavily irradiated central and 
outer regions (shaded) are in the Bose glass phase, 
while vortices in the unirradiated (white) annular 
region are in the flux liquid phase. A radial cur- 
rent drives tangential flow in the resistive flux liq- 
uid annulus, which is impeded by the "Bose-glass 
contacts" at the boundaries. The field profile ob- 
tained by solving Eq. (^ with no-slip bound- 
ary conditions |^] is spatially inhomogeneous on 
length ^, as shown in Fig. 1. One can probe 
this profile and extract ^ by placing a string of 
radial contacts at r„, for n = l,2,3,..., and mea- 
suring the voltage I4+i,n across each successive 
pair (inset of Fig. 2). If the viscosity is small 
(^ << d), the voltage decreases logarithmically 
as one moves from the inner to the outer contacts, 
as in a freely flowing uncorrelated liquid, where 
= (p///27rt)ln(r„+i/r„). When f grows, 
the onset of rigidity in the liquid becomes appar- 
ent (Fig. 2). An elastic vortex solid would ro- 
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Figure 2. The voltage drop 2TitVn+i,nl{pfI) 
across pairs of contacts (r„+i, r„), with r„ = + 
nd, for n = 0, 2, 10, Ri = d, and d = W/IO the 
contact spacing. The symbols refer to £,/d — 0.1 
(triangles), S,/d = 1 (squares) and ^/d = 2 (cir- 
cles). Solid lines are guides to the eye. Inset: top 
view of the Corbino disk with Bose glass contacts, 
with W ^ R2~ Ri. 



tate as a rigid object under the radial drive, with 



v{r) 



and V, 



{pfI/2^tRl){7 



for i?2 >> ^1- Indeed for £, > d, Vn+i,n is no 
longer monotonic with n and it exhibits a solid- 
like growth with n within a boundary layer of 
width ^. 

The resistance per unit thickness of the disk, 
PRiT,Ri,R2) = Al/(i?2,i?i)/(//27ri), is 



Pfl(r,i?i,i?2) = P/(T)^(- 



(3) 



with J^{xi,X2) a function determined by the ge- 
ometry that is obtained from the solution of 
the hydrodynamic equation. If ^ << i?i,i?2, 
viscous effect are not important and pn « 
p/ ln(i?2/-Ri) ~ 1/7. Conversely, when the 
viscous length exceeds the sample size (^ >> 
Ri, (i?2 -i?i)), then pR « P/(i?i8^C') - V??, for 
i?2 » Ri- When the effects of geometry are neg- 
ligible, resistivity measurements directly probe 
the friction, 7, while measurements in channels 
narrow compared to ^ probe the flux liquid vis- 



Table 1 

Behavior of 7 and 77 near various transitions of 
vortex matter. Here t = {T-Tq)/Tq, with Tq the 
relevant transition temperature in each case. 
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"Polymer" Glass [| 


finite 




First order freezing m 


finite 


jumps to 00 


Continuous freezing [|l^ 


finite 





cosity, 77. Experiments in artificial pinning struc- 
tures can therefore be used to infer both 7 and 
rj. The behavior of these two parameters near a 
phase transition can in fact be used to identify 
the transition, as summarized in table 1. 

Of particular interest is the case of continuous 
glass transitions, which are characterized by uni- 
versal critical behavior. Experiments of the 
type just described are especially powerful in this 
case as they can be used to map out the crit- 
ical behavior. [|j One can probe the Bose glass 
transition by lightly irradiating the liquid annu- 
lar region, and then lowering the temperature at 
constant field from T|^i"i">' < T < T^^t^cts 
T T^™"'''''. Most physical properties near 
the transition can be described via a scaling 
theory in terms of diverging correlation lengths 
perpendicular and parallel to the field direction. 



U{T) ^ IT-TbgI-"^ andC||(r) ^ \T-Tbg\~''\ 
with v\\ = 2v^, and a diverging correlation time, 
T ^ II ^ \T - Tbg\~'"'^% Scaling can then 
be used to relate physical quantities to these di- 
verging length and time scales. In particular, the 
friction coefficient 7 that determines the bulk flux 
flow resistivity Pf{T) is predicted to diverge as 
T ^ T+c as Pf ^ \T - TbgI-"^^'-^^ ■ 1 As 
shown recently by Marchetti and Nelson, when 
flux flow in confined geometries is analyzed by 
combining hydrodynamics with the the Bose glass 
scaling theory - generalized to the spatially in- 
homogeneous case - the Bose glass correlation 
length is naturally identified with the viscous 
length ^.|^ It then follows that the liquid shear 
viscosity diverges at Tbg as 77 ~ jT — TbgI^^"-^- 
Furthermore, the scaling of the finite-geometry 



4 



resistivity displayed in Eq. ^ is a general property 
of the vortex liquid near a continuous glass tran- 
sition. Hydrodynamics yields the precise form of 
the scaling function, which depends on the exper- 
imental geometry and can be found in for the 
Corbino disk. 

3. PLASTIC FLOW IN DRIVEN SOLIDS 

As shown in recent experiments by the Argonne 
group, the Corbino disk geometry can also be 
used to study the onset of plastic flow in a driven 
solid. 1^ In this case we consider an unirradiated 
disk, where the vortex array has a clear melting 
transition. Below melting the vortex solid moves 
as a rigid body, with v[r) ~ r, and the voltage 
grows as . The ^ 1/r dependence of the driving 
force yields, however, large elastic deformations of 
the medium, described by the solution of 

CeeV^u = -nQ(t>Qi. x J(r), (4) 
c 

with free boundary conditions (cge is the shear 
modulus of the lattice, assumed incompressible). 
Elastic deformations yield a finite shear stress, 
(Jrc^ir) = (no0o//4^c<)[l + 21n(i?2/i?i)i?f/r2], 
that can unbind dislocations from bound pairs. 
Assuming for simplicity that vortices are straight 
along the field direction, the energy of a 
pair of dislocations of opposite Burgers vec- 
tors, separated by a distance x, is Uq{x) — 
(ceetoo/^) [l^(^/'^o) - cos^ 6\ - 2Ect, with 9 the 
angle between x and b and E,. « cqqO^ the core 
energy per unit length of an edge dislocation. An 
applied stress pulls the two dislocations in oppo- 
site directions, [p^ Ignoring climb, and assuming 
the spatial variations of the stress field are on 
scales large compared to x ~ ao, the interaction 
energy is now U{x) = Ua{x) — aoxarcj,{r), with r 
the radial location of the center of mass of the 
pair. The applied stress lowers the barrier that 
confines the bound pair. In the simplest model, 
unbinding occurs where the location xg of the 
barrier, defined by [dU/dx'\^_^ — 0, becomes 

comparable to ao, or x_b = ceeao/ {2Trar4>ir)) ~ 
flQ. By solving for r, one obtains the crit- 
ical radius R]\i{I) where shear-induced dislo- 
cation unbinding occurs, yielding slippage of 



neighboring planes of the vortex lattice, Rm ~ 
i?i^21n(i?2/i?i)//(/o for R2» Ri, with 

Iq = 2cc66^/("'O0o) ci maximum current above 
which the entire vortex solid shear melts. The 
melting radius Rm increases with current, indi- 
cating that, since the stresses are largest near 
the axis of the disk, "shear-induced melting" oc- 
curs first in circular layers close to the axis. This 
behavior is qualitatively consistent with the ob- 
servations by the Argonne group. |^ The simple 
model described here suggests that at high fields 
the current scale Iq is independent of field. A 
more detailed calculation incorporating field and 
temperature dependence will be described else- 
where. 
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